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We discuss in detail the uniform discretization approach to the quantization of totally constrained 
theories. This approach allows to construct the continuum theory of interest as a well defined, 
controlled, limit of well behaved discrete theories. We work out several finite dimensional examples 
that exhibit behaviors expected to be of importance in the quantization of gravity. We also work 
out the case of BF theory. At the time of quantization, one can take two points of view. The 
technique can be used to define, upon taking the continuum limit, the space of physical states of the 
continuum constrained theory of interest. In particular we show in models that it agrees with the 
group averaging procedure when the latter exists. The technique can also be used to compute, at the 
discrete level, conditional probabilities and the introduction of a relational time. Upon taking the 
continuum limit one can show that one reproduces results obtained by the use of evolving constants, 
and therefore recover all physical predictions of the continuum theory. This second point of view 
can also be used as a paradigm to deal with cases where the continuum limit does not exist. There 
one would have discrete theories that at least at certain scales reproduce the semiclassical properties 
of the theory of interest. In this way the approach can be viewed as a generalization of the Dirac 
quantization procedure that can handle situations where the latter fails. 



I. INTRODUCTION 



Loop quantum gravity has emerged in recent years as one of the major candidates for a theory of quantum gravity. 
See [l| for recent reviews. The theory has a mathematically rigorous basis for its quantum kinematics [2J, which has 
also been proven to be unique up to certain assumptions [3[. It also has achieved attractive physical results in the 
context of providing a detailed microscopic picture of black hole entropy and a detailed picture of the big bang in the 
context of homogeneous quantum cosmologies. However, the problem of the dynamics of the full theory has remained 
unsettled. Related to this is the lack of a good understanding of how general relativity arises as a semiclassical limit 
of the quantum theory and the fact that the algebra of quantum constraints — though free of anomalies — does not 
mimic the algebra of the classical constraints completely. 

The reason for the elusive nature of the dynamics has a well defined but technical explanation. In the accepted and 
largely unique kinematical setup in terms of spin networks there is a natural action of the diffeomorphism group that 
is a unitary representation and therefore the spatial diffeomorphisms are represented without anomalies. However, the 
action is not weakly continuous. This means that the infinitesimal generators of the diffeomorphism group, that is, the 
associated Lie algebra, cannot be defined on the Hilbert space. In contrast, versions of the Hamiltonian constraint can 
be defined on the Hilbert space. However, since the classical Poisson algebra of constraints involves the infinitesimal 
generators of diffeomorphisms, it should not be surprising that the classical constraint algebra is not appropriately 
reflected in the quantum constraint algebra. A clear discussion of this problem is present in the introduction to the 
recent paper by Giesel and Thiemann [4J. There does not appear to exist an obvious solution to this problem in 
the sense of constructing a "better" Hamiltonian constraint operator, for obvious reasons. This obstruction has led 
several researchers to consider alternatives to the usual Dirac approach to the problem of the dynamics. One of the 
alternatives is the "master constraint" project of Thiemann and collaborators [5] motivated by an earlier proposal 
by Klauder Q. Others consider the covariant "spin foam" approach as an alternative, since one may bypass the 
construction of the canonical algebra entirely, at least in some settings. Another approach is the one furthered in this 
paper, which we will refer to as "uniform discretizations" . 

Over the last few years we have developed a new paradigm for dealing with totally constrained theories (see @, H, 
for recent reviews). The app roach has been called "consistent discretizations" and a particularly promising version 
of it recently introduced [1JJ] is called "uniform discretizations" . To describe the paradigm in a nutshell one can say 
that its relationship to totally constrained theories is the same that lattice QCD has with continuum QCD. The new 
approach consists in building a set of discrete theories and to define a continuum theory as a suitably well defined 
limit. Just like one can view lattice QCD as a tool to define continuum QCD non-perturbatively one can view the 
uniform discretization as a tool to define continuum theories that are totally parameterized, as general relativity. 
Therefore the goal is to have a discretized version of the theory that is well defined, under control and quantizable 
without difficulties and only then study how to take the continuum limit. Other approaches to the dynamics of 
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quantum gravity currently being considered also use discretizations, but the continuum limit is immediately taken 
since the discrete theories are in themselves inconsistent, i.e. they do not exist. 

The naive application of a discretization — like the ones used in lattice QCD — to gravity or other totally constrained 
field theories immediately fails since the introduction of the lattice conflicts with diffeomorphism invariance. From 
an operational viewpoint, this translates itself in the failure of constraint algebras to close. In other words, the 
discretization of a totally constrained theory generically yields a discrete theory that is inconsistent, its equations 
cannot be simultaneously solved. This is manifested by the lack of closure of constraint algebras (more precisely, 
the constraints do not structure themselves into an algebra). A typical example of such an approach is the one 
followed in numerical relativity in the Cauchy approach, where it is well known that one cannot solve simultaneously 
the evolution equations and the constraint equations. In that context the usual approach is to ignore some of the 
equations — typically the constraints — in the hope that if one starts with data that satisfies them initially, they 
will remain small upon evolution. The logic is an a-posteriori one, that is, if one carries out a simulation and the 
constraints remain small, then it is an acceptable approximation to the continuum theory. But it frequently happens 
that they do not remain small and the simulation then has to be abandoned. Much more serious is the situation if one 
is interested in quantizing. To use a set of inconsistent theories to define via a limiting process a continuum theory 
at the quantum level is clearly not logically tenable; to exist a limit one first has to have a sequence. In quantum 
mechanics one cannot just pick and choose some initial data and hope that "constraints will remain small" . Since the 
inconsistent theories do not exist, their quantization makes no sense and one cannot construct a sensible continuum 
limit. 

This is the reason why we have been studying the possibility of constructing discretizations of totally covariant 
theories that are consistent, that is, the equations of the discrete theory can be solved simultaneously. How can this be 
achieved? Largely by turning the quantities that in the continuum play the role of Lagrange multipliers into dynamical 
variables. The extra variables allow to make compatible the set of equations that were formerly incompatible. For 
instance, in general relativity written in terms of metric variables in the traditional 3 + 1 decomposition, one has twelve 
canonical variables and twelve evolution equations, plus four constraint equations. One has therefore 16 equations 
for 12 unknowns. In the continuum the equations are compatible and can be solved. When discretized however, one 
ends up with a set of incompatible equations. If we now consider the three components of the shift and the lapse as 
dynamical variables we have 16 equations with 16 dynamical variables as unknowns and simultaneous solution of the 
equations is generically possible. 

An immediate question that arises is: if the Lagrange multipliers are now determined, is this not equivalent to 
choosing a gauge? The answer is obviously negative, since the resulting theory has more degrees of freedom than the 
one we started from, it cannot simply be a gauge fixing of it. The coordinate freedom is indeed encoded in the extra 
degrees of freedom of the discrete theory. Still, since the Lagrange multipliers are determined (and in the case of a 
totally constrained theory this means the lapse is determined), how is the continuum limit to be materialized? The 
answer is that one can approximate the continuum theory by choosing initial data in the discrete theory in a careful 
way. To put it differently there exist initial data for the discrete theory such that when evolved the resulting solution 
approximates well the solutions of the continuum theory. There also exist initial data for the discrete theory that do 
not approximate the continuum at all. The continuum limit therefore is defined by choosing suitably limits within 
families of initial data. To do this in a consistent way within a canonical framework and throughout all of phase 
space is quite non-trivial. This is what the "uniform discretizations" achieve. There one discretizes the theory in a 
way that certain quantities are conserved upon evolution and therefore one has a well defined way to characterize the 
continuum limit to be taken in terms of the initial data. 

The uniform discretizations were introduced briefly in reference (Io| . The purpose of this paper is to apply the 
uniform discretization technique to a set of examples of increasing level of challenge. Most of the examples are the 
same ones as Thiemann and collaborators have considered in their master constraint program. The organization of 
the paper is as follows: in section 2 we briefly review the consistent and uniform discretization techniques, in section 3 
we consider the first example: homogeneous cosmologies; in section 4 we consider a model with Abelian constraints, in 
section 5 a model with non- Abelian constraints that structure themselves into a Lie algebra, in section 6 we consider 
a system with a non- Abelian and non-compact Lie algebra of constraints, in section 7 we explore BF theory as a first 
example of a field theory. We end with a discussion. The examples of section 4, 5, 6 have been discussed in detail 
with the "master constraint" approach Q . 
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II. CONSISTENT AND UNIFORM DISCRETIZATIONS: A BRIEF REVIEW 



A. Consistent discretizations 



We start by considering a continuum theory representing a mechanical system. Its Lagrangian will be denoted by 
L(q a , q a ), a — 1 . . . M. This setting is general enough to accommodate, for instance, totally constrained systems. In 
such case q will be the derivative of the canonical variables with respect to the evolution parameter. It is also general 
enough to include the systems that result from formulating on a discrete space-time lattice a continuum field theory, 
like general relativity. 

We discretize the evolution parameter in intervals (possibly varying upon evolution) — t n = e n and we label 
the generalized coordinates evaluated at t n as q n . We define the discretized Lagrangian as 



L(n,n + 1) = L(C<£+i) = ^L(q a ,q a ) (1) 

where 

q a = q% and q a = g " +1 - g " . (2) 

Of course, one could have chosen to discretize things in a different fashion, for instance using a different approx- 
imation for the derivative, or by choosing to write the continuum Lagrangian in terms of different variables. The 
resulting discrete theories generically will be different and will approximate the continuum theory in different ways. 
However, given a discrete theory, the treatment we outline is unique. We will take advantage of this freedom in how 
to discretize the theory to propose the "uniform discretizations" in the next section. 

The action can then be written as 

N 

S = ^L(n,n+l). (3) 

n=0 

It should be noted that in previous treatments we have sometimes written the Lagrangian in first order form, i.e. 
L = J dt (pq — H(p,q)). It should be emphasized that this is contained as a particular case in the treatment we are 
presenting here. In this case one takes both q and p to be configuration variables, and one is faced with a Lagrangian 
that involves q n ,p n and q n +\ as variables and does not depend on p n +i- The reason we frequently resort to first 
order formulations in the various concrete examples we discuss is that the Ashtekar formulation used in loop quantum 
gravity is naturally a first order one and we usually tend to frame things in a closely related way. But again, there 
is no obstruction to using either first or second order formulations with our framework, they are both contained as 
particular cases. 

If the continuum theory is invariant under reparameterizations of the evolution parameter, one can show that the 
information about the intervals e„ may be absorbed in the Lagrange multipliers. In the case of standard mechanical 
systems it is simpler to use an invariant interval e„ = e. 

The Lagrange equations of motion are obtained by requiring the action to be stationary under variations of the 
configuration variables q a fixed at the endpoints of the evolution interval n = 0, n = N + 1, 

dL{n,n + \) dL(n-l,n) = 

dq- + dtf* ■ { > 

We introduce the following definition of the canonically conjugate momenta of the configuration variables, 

dL(n,n+l) 

Pn+1 ~ (5) 
a = dL(n - l,n) = dL(n,n+l) 
Pn " dql dq« { > 

Where we have used Eq. (|4j) . The equations (|5|) and (|6j) define a canonical transformation for the variables qn,Pn 
to q n +i,p n +i with a the type 1 generating function iq = — L{qf ll <7° +1 ). Notice that the evolution scheme is implicit, 
one can use the bottom equation (since we are in the non-singular case) to give an expression for q n +i in terms of 
q n ,p n , which in turn can be substituted in the top equation to get an equation for p n +i purely in terms of q n ,p n - 
It should be noted that there are several other possible choices, when going from the set of equations (|5l6p to an 
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explicit evolution scheme (see Di Bartolo et al. [Tl| for further details.) 
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The canonical transformation we introduced becomes singular as an evolution scheme if 



d 2 L(n,n+ 1) 



vanishes. 



If the rank of the matrix of second partial derivatives is K the system will have 2(M — K) constraints of the form, 

<MM) - o (7) 

= 0. (8) 

And these constraints need to be enforced during evolution, which may lead to new constraints. We refer the reader 
for the detailed Dirac analysis to Di Bartolo et al. fill ]. 

To clarify ideas, let us consider an example. The model consists of a parameterized free particle in a two dimensional 
space-time under the influence of a linear potential. The discrete Lagrangian is given by, 

L n = L(q^,TT^, N n ,q^ +I ,n^ +Il N n+1 ) (9) 

= <(Qn+i - O - NnK + \K) 2 + aqi\. 

We have chosen a first order formulation for the particle (otherwise there are no constraints and the example is trivial). 
However, this Lagrangian is of the type we considered in this presentation, one simply needs to consider all variables, 
q a , 7r a , N as configuration variables. The system is clearly singular since the 7r's and N only appear at level n (or in 
the continuum Lagrangian, their time derivatives are absent). When considered as a Type I generating function, the 
above Lagrangian leads to the equations 

K,n + 1 = ^=0, (10) 
un n+l 

" Ln ' (11) 



J q, n+1 



dq a n+1 



PN -^ = Mtr^ (12) 



and 



= = ~ €) + <N n 51 + N n 5 a , (13) 

Pin = ~^=K+^aN n , (14) 

Pn,u = -Qjy- = tt„ + -(tt„) +aq n . (15) 

The constraints (|10I12I14I15P can be imposed strongly to eliminate the 7r's and the iV's and obtain an explicit evolution 
scheme for the g's and the p q 's, 





- n° 


Cn+1 
"Pj, n+l ' 


(16) 


Qn 


— In+l 


Cn+1 

a 


(17) 


i r g, n 


- n° 


-l) 


(18) 


Pq. n 


= Pq, n-i 


, C n+ i 
-i+i i 

Pq, n+l 


(19) 



and the Lagrange multipliers get determined to be, 

N n = (20) 

CXPq, n+l 

where C n+ i = jp q n+1 + (pj n+1 ) 2 /2 + aq^ l+1 . The evolution scheme runs backwards, one can construct a scheme that 
runs forward by solving for N and ir at instant n when imposing the constraints strongly. The two methods yield 
evolution schemes of different functional form since one propagates "forward" in time and the other "backward" . The 
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inequivalence in the functional form stems from the fact that the discretization of the time derivatives chosen in the 
Lagrangian is not centered. It should be emphasized that if one starts from given initial data and propagates forwards 
with the first system of equations and then backwards using the second, one will return to the same initial data. 

So we see in the example how the mechanism works. It yields evolution equations that usually are implicit as 
evolution schemes. The equations are consistent. The Lagrange multipliers get determined by the scheme and 
there are no constraints left on the canonical variables. The evolution is implemented by a (non-singular) canonical 
transformation. The number of degrees of freedom is larger than those in the continuum. There will exist different 
sets of initial data that lead to different solutions for the discrete theory but nevertheless will just correspond to 
different discrete approximations and parameterizations of a single evolution of the continuum theory. 

The example also exhibits some of the problems the framework may face. Since one is generically left with implicit 
non-linear systems of algebraic equations to be solved, there is no guarantee that the solutions will be real. Neither 
is it guaranteed that the Lagrange multipliers will be bounded in value (this is important since in totally constrained 
theories the lapse controls the evolution step and one wishes it to be small in the continuum limit). Although we 
explored some simple examples where this does not pose a problem in several publications, it is clear that generically 
it could be a serious obstacle. This obstacle can be related to the intuitive fear that people manifested about the 
formalism being a "gauge fixing" . Gauge fixed quantizations are not favored since it is known that one cannot fix 
gauges globally. Similarly here, the framework cannot generically guarantee that one approximates the continuum 
globally in a solution. We will now argue that some simple, yet carefully chosen, discretizations actually achieve this 
goal. This is in principle achieved in a generic fashion, i.e. we can find such discretizations for any given totally 
constrained system. 



B. Uniform discretizations: the key idea 

As we stressed before, there is significant freedom in how one chooses to discretize the action. We will now use 
this freedom to our advantage, as was first proposed in Consider a totally constrained system with N first class 
constraints fa and a configuration space with 2P variables (q,p) (we omit indices on the variables for simplicity). We 
will assume that one has chosen the discretization of the action in such a way that the evolution equations for a given 
dynamical variable A of the theory can be written as, 

A n+1 = e^(A n )=A n + {A n ,H} + ^{{A n ,H} 7 H} + --- (21) 

As is obvious, H is a constant of the motion of the discrete evolution. The quantity H is defined in the follow- 
ing way: Consider a smooth function of N variables f{xi, . . . , xn) such that the following three conditions are 

satisfied: a) f(xi, . . . ,xn) = <^=> xi = OVi and otherwise / > 0; b) J^(0, ...,0) = 0; c) det d ®J x . ^ 

\fx and d) f(<pi(q,p),...fa\j(q,p)) is defined for all (q,p) in the complete phase space. Given this we define 

H(q,p) = f(4>i(q,p), ■ ■ ■ (j) N {q,p))- 

A particularly simple example is H{q,p) — 1/2^3- j fa(q,p) 2 , a choice that has interesting parallels with the 
"master constraint" of the "Phoenix project" [5| as we shall discuss later. 

An important point is that if we choose initial data such that H < e then fa remain bounded throughout the 
evolution and will tend to zero in the limit e-> 0. Let us see that in this limit one recovers the evolution equations 
given by the total Hamiltonian Ht in the constrained continuum theory. Let H as in the simple example above 
and take its initial value to be H = 5 2 /2. We define Aj = fa/5, and therefore J2iLi = 1- The evolution of the 
dynamical variable q is given by, 

N 

qn+i = q-n + J2{^fa}KS + 0(S 2 ) (22) 

i=l 

and if we define q = lim^o (in+i ~ Qn)/S, where we have identified the "time evolution" step with the initial data 
choice for 5, one then has, 

N 

q = Y,{q,fa}Ai, (23) 
i=i 

and similarly for other dynamical variables. The specific values of the multipliers Xi depend on the initial values of 
the constraints fa . Notice that taking the continuum limit requires that the Lagrange multipliers be determined as is 
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usual in the consistent discretization approach, but are well defined bounded real functions of phase space, bypassing 
an important objection to the original approach. 

Notice that the procedure deals only with first class constraints. In case of a system with second class constraints, 
it will be needed to use Dirac brackets and work with the equivalent first class system. 

An important comment is needed here about the difference in the case of a mechanical system and a field theory. 
Although it is true that field theories formulated on the lattice are mechanical systems, generically the constraints of 
field theories formulated on the lattice fail to be first class, even if they were so in the continuum. If the constraints of 
the discrete theory are not first class the above procedure fails to reproduce the equations of the continuum. Here one 
has two choices. One would be to first take the continuum limit spatially keeping time discrete. In that case it may 
happen that the resulting constraints of the theory with continuum space and discrete time are first class. Then the 
above proof goes through. The alternative is to consider in the discrete evolution equation the Dirac bracket of 
the discrete theory. Then one can show that the continuum limit is correctly achieved, including in the Hamiltonian 
both first and second class constraints. 

One could recast the current proposal in terms of the original approach to consistent discretizations. There one 
started from an action and noted that the Lagrangian could be viewed as a the generating function of a canonical 
transformation between instants n and n + 1. To be concrete, let us analyze a system with N Abelian constraints. 
We introduce the Lagrangian 



where L is a type 1 generating function of a canonical transformation between canonical variables q n , p n and q n +\ , p n +i , 
S is Hamilton's principal function for a given set of Lagrange multipliers Ai , . . . , A at (they are evaluated at instant n, we 
omit the subscript for simplicity), g is such that g(0) = and the mappings Ai — > Jj- and Xi — > Jj- are inverse where 
/ is the function used to define the Hamiltonian. The generating function yields the canonical momenta in the usual 
way Pn+i = 9L/dq n +ii Pn = —dL/dq n . One also has that dL/dXi = and this determines the Lagrange multipliers, 
A; = hi((p), where hi is the inverse function of the mapping defined by Ai — > Jj-. This evolution corresponds to a 
Hamiltonian H — f(<f>i, ■ ■ . , 4>n), with dif = hi. In particular, the simplest case is when g = X)t=i x i/% an d then 
H = Yli=i The generating function L allows to determine the discrete evolution that preserves exactly the 

value of the constraints of the continuum theory and recovers the continuum limit when all <jn — ► in the initial data. 
We therefore see that the approach proposed here is a particular case of the consistent discretizations. We have just 
chosen to discretize things in a way that the Hamiltonian is simple — rather than the action — and this guarantees a 
good continuum limit. 

The constants of the motion of the discrete theory are quantities that have vanishing Poisson bracket with the 
Hamiltonian, {Of , H} — and in the continuum limit Ho — > reproduce, as functions of phase space the "perennials" 
of the continuum theory Of — lim.ff D _>o Of . This can be immediately seen from the fact that the discrete equations 
reproduce the continuum equations for any dynamical variable in the continuum limit. Conversely, for every perennial 
of the continuum theory there exists a constant of the motion (in general many constants) of the discrete theory that 
reduce to the given perennial in the continuum limit. We have therefore shown that uniform discretizations recover the 
constraints and the perennials of the continuum theory and therefore provide a good starting point for a quantization 
of the continuum theory. 



We now turn our attention to the quantum theory. We will introduce a Heisenberg quantization for the discrete 
theory (this is more natural given that one has an explicit evolution). To quantize the theory we follow several steps. 
We start with the classical discrete system constructed as in the previous section, we then have the canonical variables 
at level n + 1 in terms of the variables at level n, q n +i = Qn+i(qn,Pn), Pn+i — Pn+i(qn,Pn)- 

We then define the kinematical space of states of the quantum theory, Hk, as the space of functions of N real 
variables ^(q) that are square integrable. In this space we define operators Q and P as usual. To construct the 
operators at other time levels (in the Heisenberg Picture) we introduce a unitary operator U that we will define later 
that gives, 



L(q n , q n +i, Ai, . . . , Xn) = S(q n , q n +i,\i, . . . , Xn) + g(Xi, . . . , Ajv) 



(24) 



C. Quantization 



Pn = 



u- 1 P n - 1 u = u- n p u n . 



(25) 



When the evolution is determined by a discrete Hamiltonian H, as is the case in the uniform discretizations, the 
evolution operator is given by U — e~ lB / h . Notice that U may also be determined by requiring that the fundamental 
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operators satisfy an operatorial version of the evolution equations, 

QnU-UQ n+1 (Q n ,P n ) =0, P n U -UP n+1 {Q n ,P n ) = 0, (26) 

and this provides a consistency criterion for the construction of U. 

At a classical level H = if and only if the constraints fa = 0. There exists a natural definition of the physical space 
of the continuum theory that does not require that we refer to the constraint. Since we know that U = exp(—iH/h), 
a necessary condition satisfied by the states of the physical space of the continuum theory, tp S H p h ys is given by 
Uip — tp. More precisely the states ip of H p h ys should belong to the dual of a space <I> of functions sufficiently regular 
on Ji k . That is, the states ip G W p h ys satisfy 

ip*U^ipdq = [ ip*tpdq, (27) 



where tp S $. This condition characterizes the quantum physical space of a constrained continuum theory without 
needing to implement the constraints as quantum operators by using the discretization technique. 

The unitary operators of the discrete theory allow to construct the "projectors" onto the physical space of the 
continuum theory, which is one of the main goals of any quantization procedure based on Dirac's ideas. It should be 
noted that these are really generalized projectors in the sense that they project to a set of functions that belong in 
the dual of a subspace of sufficiently well behaved functions of Hk- All of this is achieved without having to define 
the quantum constraint. To construct the "projectors" one can compute, 

P= lim C M U M . (28) 

M— >oo 

If such a limit exists for some Cm such that ]im.M^>oo{C m+i / C m) = 1 then UP = P, and we have that UPip = Pip, 

yip s Hk- 

The limit exists in several examples in which H has a continuum spectrum, as we shall see. If the spectrum is discrete 
with eigenvalues ej and it contains a vanishing eigenvalue ei then a projector is trivially defined as |ej >< ej |. A 
constructive procedure leading to a general definition of the projector in terms of the discrete evolution operator U 
valid for any spectrum, continuum or discrete, is given by: 

Int(rM) f 

lim y 7—777 77 - (29) 
M -oo ^ Int(rM) — M K ' 

n=M v ; 

where r is a real number grater than one and Int(rM) is the integer part of rM . If U has a continuum spectrum 
this definition is a trivial consequence of the previous one. In the case of a discrete spectrum one can check that the 
definition works recalling the definition of the Kronecker delta in terms of a Fourier series. Notice that the definition 
of physical space that Thiemann introduces in the "Phoenix project" is equivalent to the choice we make if one 
is considering the Hamiltonian that is quadratic in the constraints. Furthermore, given two states of 7i p h ys , V'ph, 
(f> p h, where ip p h = Pip, and </> p h = P<f>, the physical inner product is defined by < ip p h\<pph >= f dq<fi(q)* Pip(q) and a 
physical inner product is determined by the projector constructed from the discrete theory. 

We now illustrate the technique with a rather general example. We consider a generic mechanical system with 
a finite dimensional phase space with one constraint = 0. We will show that the projector constructed with our 
technique reproduces the one constructed with group averaging techniques [121 ] . That is, 

P= lim .[JK e -iMt* = r oo) 

Af^oo V TT J_ OQ 27T 

To make contact with the group averaging case we need to assume that is a self-adjoint operator in the kinematical 
phase space with an eigenbasis given by <fi\a >= (p{a)\a > and 1 = J \a >< a\da. The proof of the equivalence is, 




P = P \a><a\da= / lim W e~ %M,p \a >< a\ (31) 



and noting that limAf^oo J ^-e lMx = S(x) = ^e ltix the proof is complete. For the proof we assumed a 
quadratic form of the Hamiltonian, but can actually be extended to Hamiltonians of the general form we discussed 
above, computing the integral by steepest descents. The proof can also be extended to systems with TV Abelian 
constraints by noting that limM^ 1 ^)^-^-* = 5{x) = J'^ d^ N e 1 ^ /(2w) N . 

A point to be noted is that if the Hamiltonian is such that zero is not an eigenvalue, the above expression for the 
projector vanishes. Therefore the technique cannot be applied. We will confront a situation like this in section VI. 
There we will argue that in such models one should take a different point of view than defining the physical Hilbert 
space in the traditional way. 
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III. EXAMPLE 1: HOMOGENEOUS ISOTROPIC COSMOLOGY 

We will follow the same procedure as in ordinary loop quantum cosmology (see [l3T | for a review). We start by 
considering general relativity formulated in terms of Ashtekar's variables. We then assume that we have a homogeneous 
and isotropic universe with flat spatial sections. In such case one can choose the triad and Ashtekar connection as, 

cS { v8 a 
A IW = $fc> E ? = ( 32 ) 

with {c,p} = (87r/3)G where we have taken the Immirzi parameter as unity. As usual one proceeds in analogy with 

1 /3 

the full theory and defines the holonomy along an "edge" of "length" p,V , where Vb is the volume an elementary 
cell of universe, 

(c) = cos(^c/2) + 2r, sin(/xc/2) . (33) 

The loop quantum cosmology approach, roughly speaking, can be characterized as attempting to produce for the 
homogeneous context expressions that mimic those that are found in the loop representation for the full theory. In 
particular, one takes the usual Hamiltonian constraint of the classical theory and represents the quantities involving 
connections using holonomies and the expressions involving the triads using the Poisson bracket of the volume with 
the connection, as first noted by Thiemann in the full theory 

The resulting expression for the Hamiltonian is, written as 

°^ = £ ^ (^ 0) ^° ) ^ 0) " 1 4 M0) K 0) ~\ V}) + O(cVo), (34) 

ijk 

where V = \p\ 3 ^ 2 - In the continuum limit (fio — ► 0) it can be shown that the expression reproduces the Hamiltonian 
constraint of a homogeneous isotropic vacuum minisuperspace. 

It turns out to be more convenient to evaluate explicitly the Poisson bracket and replace the holonomies by the 
expression above to get, 



C <£) = -V6 Vl6 f lpl sin 2 ( Mo c) . (35) 
Mo 

Where we have made the redefinition -^§qP — > p so that {c,p} = 1, we will use this redefinition of p in what follows. 
This Hamiltonian constraint can be obtained from a Lagrangian for gravity coupled to a scalar field of the form, 



1 M 
L =pc + p,p<p+ — 



|p|sin 2 (^ c) P\ 



Mo bl 3/2 



(36) 



M is the Lagrange multiplier (proportional to the lapse) with the factor V6Vl67rG absorbed in it. For simplicity we 
have also rescaled the scalar field and its conjugate so the gravitational coupling 3/(16ttG) is absorbed. 

We immediately recognize in the classical limit [io —> the usual Lagrangian for the Friedmann-Robertson- Walker 
model coupled to a scalar field. 

Given that we are dealing with a mechanical system with one constraint, one can immediately apply the technique 
used in the example of the last section to show that one recovers the same physical space as in the group averaging 
approach. We therefore make contact exactly with traditional loop quantum cosmology, at the level of physical space. 

It would be interesting to pursue this example further by not taking the continuum limit. There one could introduce 
a relational time and study loss of unitarity and other effects that will arise, but we will postpone a detailed discussion 
of cosmology to a subsequent paper. 



IV. EXAMPLE 2: A MODEL WITH ABELIAN CONSTRAINTS 

We consider a mechanical system with configuration manifold R N , coordinatized by q l , i — 1, . . . , N and M < N 
commuting constraints 



C l =p l i = l,...,M 



(37) 
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where the p l, s are the conjugated momenta to the q l, s. 

All phase space functions which do not depend on the first M q l, s are Dirac observables, i.e. functions which 
commute with the constraints (on the constraint hypersurface) . 

A Dirac observable which depends on p', i = 1, . . . ,M is equivalent to the Dirac observable obtained from the 
first one by setting p l = (since these two observables will coincide on the constraint hypersurface). Therefore it is 
sufficient to consider observables which are independent of the first M configuration observables and of the first M 
conjugate momenta. A canonical choice for an observable algebra is the one generated by q l ,p l , i = (M + 1), . . . , N. 

Let us now consider the uniform discretization. For that we construct H = | YLiLi (p 1 ) 2 - The constant k is to 
make the Hamiltonian a dynamical variable with units of action. The discrete evolution equations are, 

M 

= ?n+fc£{e^R, (38) 

6=1 

= + fcp°9 QiM , (there is no summation in a) (39) 

Pn+l = Pn> (40) 

where O a ,M = 1 if a < M and zero otherwise. 

To take the classical continuum limit, we assume that H = S 2 /2 with S — > 0. If we define 



we have from the evolution equations that, 



■a = Um gn+1 gn ^ 
8^0 S 



q a = A a 9 aM , (42) 
P a - 0, (43) 



where A a = lim^^ kp a /S. 

To quantize the system, we will start with an auxiliary Hilbert space L2{R N ) on which the operators <f act as 
multiplication operators and the momenta as derivatives, 



and construct a unitary operator, 



in terms of which one has 



= (44) 
pty(«) = -iKdrt(q), (45) 



C^ = exp^-iA^;^, (46) 



C+i = (47) 
M k 

= C + E^[(P«) 2 '^] (48) 

6=1 

= q% + kp a n Q aM , (49) 

so we see that we have recovered the classical equations as operatorial relations in the Heiscnberg representation. 

We now wish to define the physical space of states of the continuum theory. To do this we will use the projector 
technique, 

P= lim C R U R . (50) 

R—y oo 

with Cr — (|ff ) M ^ 2 - Acting on the identity written as I — J dp 1 ■ ■ ■ dp N \p 1 , . . . ,p N >< p 1 , . . . ,p N \ we have, 

pi= 1'ip' - ->ip s ib" o,.oxpi -'--y (//')-' I \i>'....,p N ><p 1 ,..., P N \ 



1 



J dp^.. dp N R lim C R exp (-i^fy*)^ \p\ 



dp M +i ...dp N \0,... ,p M +\ ...,p N ><0,... ,p M +\ ...,p N \. (51) 
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We therefore see that the physical space is given by the square integrable functions depending on the variables 
M + 1 to N. The physical inner product is therefore immediately induced by the kincmatical inner product on the 
physical space of states. 

We would like now to turn our attention to a quantization where we assume the continuum limit is not taken, i.e. we 
end up with a discrete quantum theory. Although this is not necessary in this example, it would be in examples where 
the continuum limit is problematic. It is therefore good to see how the discrete quantum theory can be viewed as a 
quantizations of the continuum theory, at least at certain scales. Since we have a conserved Hamiltonian if we choose 
an initial quantum state that has a small expectation value for the Hamiltonian, it will keep its small expectation 
value upon the unitary evolution defined by U. We will take advantage of working in the discrete theory to be able 
to measure variables that would not be observable in the continuum theory and in terms of them define a relational 
evolution that will approximate well the dynamics of the continuum theory. 

We start by recalling one way in which the dynamics of models like this can be handled in the continuum via the 
use of "evolving constants" (la |. Given any function of the dynamical variables of the theory f{q 1 , . . . , q N ,p x , . . . ,p N ) 
one can define an "evolving constant" . This is a function that depends on a given number of real parameters 
/(T 1 , . . . , T M , q M+1 , . . . , q N ,p M+1 , . . . ,p ) and that coincides with the quantity of interest on the constraint surface 
when the partial observables q 1 , . . . , q M take the values T 1 , . . . , T M . The "evolving constant" is a physical observable 
of the theory, that is, it has vanishing Poisson brackets with the constraints. To determine the dynamics of the 
theory is to determine what is the value of the quantity of interest / for a given value of the parameters T 1 , . . . ,T M . 
A potential drawback of this approach is that it requires to treat some of the variables of the problem as classical 
variables, something that may not appear very natural in certain regimes where all variables should exhibit quantum 
mechanical behaviors. 

We now go back to the discrete theory. Since the theory has no constraints we can work in the kinematical space and 
all variables will be quantum mechanical in nature. We wish to describe the evolution relationally, defining conditional 
probabilities (see for instance fl6jV The physical variable we will consider is F = q 1 p N . This operator coincides with 
the angular momentum Lin — q 1 ^ — q N p x evaluated on the constraint surface. We now need to choose a "clock". 
Here we need to be careful. If we choose q 1 , . . . , q M as "clock" variables the problem that may arise is that if one 
makes a very precise measurement of one of these variables would imply a large change in the expectation value of 
the Hamiltonian that would lead to lose the regime close to the continuum. This is understandable, in a theory with 
fundamental discreteness there has to be a limit as to how accurate a measurement one can make. For instance, this 
would be akin to claiming that in a discrete theory of quantum gravity one cannot measure lengths shorter than the 
Planck length. It should be recalled that this is also a rather artificial limitation. If one is dealing with a realistic 
theory that includes measuring devices in it, since the evolution is unitary, one is guaranteed to remain in the regime 
one chooses irrespective of what measurements are made with the devices the theory incorporates. One can only 
break the regime if one introduces measuring devices that are external to the theory. In the particular example we 
are considering the theory is too simple to encompass measuring devices within it. We are therefore forced to consider 
external measurement devices and therefore we need to limit their interactions with the model to respect the choice 
of initial conditions within the regime of being close to the continuum. 

Taking the last point into account, the projector for a clock reading a value (q 1 , . . . q M ), where we understand that 
each variable takes a value within an interval Aq centered at each q a , is 



P 




(52) 



and the projector for a measurement of F to be in an interval (/i, f 2 ) is, 



Pf 




. . . , q M ,p M+ \ ...,p N ><q\..., q M ,p M+ \ . . . ,p N \9(f 1 /q\p N J./q 1 ) (53) 



where 9(a, b, c) — 1 if a < b < c and zero otherwise. 

The projector for a simultaneous measurement of the clock and F therefore is, 



P q P F = I dq l ■ ■ ■ dq M [ • • • dp N \q\ . . . , q M ,p M+ \ . . . ,p N X q\ . . . , q M \p 




,M+1 



...,p N \9(f 1 /q\p N J 2 /q 1 ). 

(54) 



Let us consider a state at instant n = of the form, 




(55) 
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therefore, upon evolution, 

\i, n> =U n \^>= J d p 1 ---dp N c p , ! ... !pN e X p(-^J2^ 2 ] \P\---,P N >- (56) 

The conditional probability is given by V(F\q) = (J2 n < ^ n \P q PF\^n >)/E n < ^n\P q \^n >) that is, the ratio 
of the probability that the clock measures q and the system under study takes the value F simultaneously and the 
probability that the clock measures q [l6[ . Explicitly, 

[ W 2 n Ja do 1 '" dq M J dp M+1 ■ ■ ■ dp N \D n (q 1 , . . .,q M ,p M+1 , . . . ,p N )\ 2 [b '> 

where we have defined 

M 



D n (q\...,q M ,p M +\...,p N )^ J dpi---dp M c pl _ pN f[ <q J\pi > e X p(- l ^r 



(58) 



If we assume that the initial state at n = takes a tensorial product form (which means the "clock" is not interacting 
with the "system under study") c p i^ p n — a p i p Mb p M+i „n we have 

/ A rfg 1 - dq M J ■ ■ ■ dp N 9(f 1 /q\p N , f 2 /q 1 )\b p M +1 _ pN \*h(q\ . . . , q M ) 

J A dq 1 ---dq M J dp M + 1 ---dp N \b p M+^,„ tp N\ 2 h(q 1 ,...,q M ) [ ' 

for the conditional probability. Here we have defined 



%\...,^)-E 



f dp 1 ■ ■■dp M a p x t ... iP M Y[ < qJ\ P 3 > exp (~$V') ; 



(60) 



Since we are interested in comparing the results of our approach with that of the evolving constant approach in the 



continuum, we will now take the continuum limit. It is given by, a p i ... tP M — dip 1 ) ■ • ■ S(p M ) the ^.(q 1 , . . . , q M ) become 



constants and the probability is 

v(m s = L I d P M+1 ■ ■ ■ dp N 0(f 1 /q\p N ,h/g 1 )\b(p M +\ ■ ■ ■ ,p jV )| 2 
{W ■■■dp N \b{p M +\...,p N )\ 2 1 ' 

This completes the calculation in the relational picture. 

Let us now revert, for comparison, to the "evolving constant" picture. To do this we need to construct a state of 
the physical space of the continuum theory such that it yields the same expectation value for the observables that 
\ipQ>. It is given by a density matrix, 

ft>hys = Tr|gi...g M |V>o >< (62) 

The physical quantity to measure is F = q lr p N (where q 1 is viewed as a c-number and therefore this quantity is an 
observable). The probability of getting a value in the interval (/i, f 2 ) is 

v ^ Tr (J A dqlfdpl---dp N \p\...,p N ><p\... ,p N \9(f 1 /q\p N , h/q 1 )^) 

ATr (pphys) 

= f A dqifdp^ 1 ■ ■ ■ dp N e(f 1 /q\p N , f 2 /q 1 )\b( P M +\ . . ■ ,p jV )| 2 

Afdp M + 1 ---dp N \b(p M + 1 ,..., P N )\ 2 [ ' 

and in the case in which the state is a tensor product of the form c p i ) ,,, <p n = dpi,... pMb p M+i t ,„ tP N this probability 
indeed coincides with (|6"T|) . 

We have therefore seen how one can work in the discrete theory to introduce a relational notion of time in the 
quantization of the discrete theory and how the physical predictions of this picture coincide with those of the "evolving 
constants" picture in the continuum theory if one takes the continuum limit. The agreement is qualified in the sense 
that we cannot carry out measurements of arbitrary precision of "classical" parameter of the evolving constant. This 
is to be expected in a treatment where all variables are treated quantum mechanically in an equal footing without 
singling out one to behave classically. 
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V. EXAMPLE 3: A SYSTEM WITH COMPACT GAUGE GROUP AND CONTINUUM LIMIT 



Here we will consider a system with non-Abelian constraints that form a Lie algebra. The model is given by 
the same phase space as in the previous section restricted to three configuration space dimensions. We consider a 
mechanical system with configuration manifold R 3 , coordinatized by q l ,i — 1,2,3 and 3 non-commuting constraints 

Qi = L i = e iifcgjpfc i = l, 2, 3 , (65) 

where the p"s are the momenta conjugate to the q"s and we assume Einstein's summation convention on repeated 
indices. Notice that the three constraints are not independent. Vanishing angular momentum is equivalent to requiring 
q % — Xp l with A an arbitrary number, which implies two conditions between the phase space variables. The constraint 
surface is four dimensional and the system will have two independent observables. In our approach we could choose 
to construct the Hamiltonian starting from two independent constraints or one could choose to use a more symmetric, 
yet redundant, form of the constraints. We will choose the latter. We introduce a canonical evolution through a 
Hamiltonian, 

»=1 

where k is a constant with units of action. The discrete evolution equations are, 

i i q LP m L n L (an\ 

q n+ i = q n + e i]k— \-eijkej mn q — p h • • ■ (.07 J 

p k L^ L n L k 

Pn+1 = P l „ + ^jk—^+£ l jk£ 3 mnP m —£2- + --- (68) 

The continuum limit is obtained setting H = S 2 /2 and defining A 4 = L 1 / (Syk) we have, taking the limit as explained 
before, 



q k X j 



q* = eijk^-fr, (69) 



p k \ j 

We note that the components U are constants of the motion and therefore X are constant. There exist three 
further constants q-q — q l q l , p-p =p t p i , q-p = q l p l - These are not independent since L 2 — {{q- q)(jp-p) — (q-p) 2 ) = 0. 
One therefore has five independent constants of the motion of the discrete theory. In the continuum limit the L"s 
vanish and one is left with two independent constants of motion, for instance q ■ q and q-p. In the continuum theory, 
the trajectories are arbitrary trajectories on a sphere. In the discrete theory, when one takes the continuum limit 
one obtains trajectories that correspond to arbitrary circumferences on the sphere, since the A"s are constant. The 
constraint surface is therefore completely covered, but not all orbits of the continuum theory. This however, does not 
cause problems since we can recover all physical information of relevance with the trajectories we get. 

To quantize the model, we construct the same kinematical space as in the previous example (restricted to a three 
dimensional configuration space), and immediately construct the unitary operator, 



U = ^P -te . (71) 




which recovers as operatorial relations the classical discrete evolution equations (up to terms of order h) so the 
correspondence principle is satisfied. 

To compute the projector we use the formula (|29|) and, in terms of the a basis labeled by the radial and angular 
momentum eigenvalues, 

P\n,£,m >= 8io\n,£,m >= <^o|ra,0,0 >, (72) 

which can be rewritten as, 

oo 

P = K°,0 >< 0,0, n|. (73) 

n=0 
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The continuum limit is therefore trivially achieved. The physical space is the space of zero angular momentum and 
there the constraints of the continuum theory are automatically satisfied. We therefore recover the traditional Dirac 
quantization. The physical space is given by the square integrable functions depending on the radial variable \q\. The 
physical inner product is therefore immediately induced by the kinematical inner product on the physical space of 
states. 

We will not work out in detail the quantization in which the continuum limit is not taken, since it runs along similar 
lines as what we discussed in the previous section. The only thing to notice is that the current model would have a 
"time variable" consisting of two angles in the sphere. Again, one has restrictions on the accuracy with which one 
can measure angles so the approximation of the continuum theory is not lost. 



VI. EXAMPLE 4: A SYSTEM WITH NON-COMPACT GAUGE GROUP 



Here we will consider a system with non-Abelian constraints that form a Lie algebra that is non-compact. More 
specifically, the model is associated with the SO(2, 1) group. It has the same phase space as in the previous section, 
except that the metric will be given by diag(— 1, 1, 1). The configuration variables are q % with i = 0,1,2 and the 
conjugate momenta are pi. The constraints are 

C t =L l = e lJ k q3p k i = 0,l,2 , (74) 

we assume Einstein's summation convention on repeated indices. Here one has to be careful about upper and lower 
indices since the metric is non-trivial. As in the previous example we notice that the three constraints are not 
independent. The constraint surface is four dimensional and the system will have two independent observables. In 
our approach we could choose to construct the Hamiltonian starting from the two independent constraints or one 
could choose to use a more symmetric, yet redundant, form of the constraints. We will choose the latter. We introduce 
a canonical evolution through a Hamiltonian, 

H= LiLi + 2Ll , (75) 
where L is a constant with units of action. The discrete evolution equations are, 

Qn+i = In + ^jk— 2e Qk — — + ... (76) 

v k U n k L a 

Pn+l = K + t l 3 k— - ^Ok— + ■■■ (77) 

The continuum limit is obtained setting H = S 2 /2. Defining A 1 — L l /(dVk) for i = 1,2 and A = —L°/(5Vk) we 
have, taking the limit as explained before, 

<r = e ^^7p (78) 

V 1 = (79) 

We note that although Lq is a constant of the motion, the other components are not. However L 2 + L 2 is a constant, 
so L\ and Li rotate around Lq throughout the evolution. There exist three further constants q- q = q l qt, p ■ p = p l Pi 
q ■ p = q l Pi- These are not independent since L ■ L — ({q ■ q)(p ■ p) — (q ■ p) 2 ) = 0. One therefore has four independent 
constants of the motion of the discrete theory. In the continuum limit two of these vanish and one is left with two 
constants of motion, for instance q ■ q and q ■ p. In the continuum theory , the trajectories are arbitrary trajectories 
on a two hyperboloids, one space-like and one time-like. In the discrete theory, when one takes the continuum limit 
one obtains trajectories that correspond to particular choices of the Lagrange multipliers, depending on the initial 
conditions chosen. 

The quantization of this model (and similar ones) is known to be problematic (l7l . fTl . fig ]. The core of the 
problem is that if one promotes the constraints to operators on the usual Hilbert space, they do not have a vanishing 
eigenvalue. This can happen, but usually the resolution is to extend the Hilbert space by including an improper basis 
of eigenvectors. This can be done in this case, but one finds that the spectrum again does not contain zero. More 
specifically, the continuum spectrum has eigenvalues larger than H 2 /4. One can find eigenvectors with eigenvalues 
smaller than H 2 /4, but they do not arise as limit of functions of the Hilbert space, i.e. they cannot form part of an 
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improper basis of the Hilbert space. One can adopt the point of view that nevertheless the eigenvectors with zero 
eigenvalue are the "physical space" of the theory, essentially abandoning the idea that the physical space arises as 
a suitable limit of the kinematical space. This was advocated by [13, Esj. From the point of view of our technique 
this is not satisfactory, since we wish to construct the physical space of states starting from the quantum kinematical 
space, taking a limit. As we mentioned before, lacking the zero eigenvalue in the spectrum of the Hamiltonian yields 
our prescription for the projector on the physical space useless. With this in mind, the most satisfactory solution is 
the one chosen in [l8|], where one chooses as physical space the eigenvectors that have Ti 2 /4 as eigenvalue. Another 
attractive possibility in the discrete approach is not to take the continuum limit and retain a level of fundamental 
discreteness. This is very natural in a model where it is difficult to achieve a vanishing eigenvalue for the constraints, 
a natural minimum existing for their eigenvalues. As we argued before, such models can approximate the semiclassical 
physics of the theory of interest with some restrictions on the type of states considered. There can therefore be viewed 
as the best thing one can do in terms of having an underlying quantum theory for the model that approximates the 
classical physics of interest. 



VII. GRAVITY IN THREE DIMENSIONS 

In this section we will follow closely the construction of Noui and Perez [2(| • We will therefore not repeat all the 
details here. The difference will be that they construct the physical space of states using group averaging and we will 
use the projector using the technique we discussed earlier in this paper. 

We start with three dimensional gravity in first order formalism. The space-time M. is a three dimensional 
oriented smooth manifold and the action is simply given by 

S[e,u] = I Tr[e A F(u)] (80) 
Jm 

where e is the triad, i.e. a Lie algebra valued 1-form, F(u) is the curvature of the three dimensional connection ui 
and Tr denotes a Killing form on g. For simplicity we will concentrate on Riemannian gravity so the previous fields 
should be thought as defined on SU{2) principal bundle over M. 

The phase space is parameterized by the pull back to S of u and e. In local coordinates we can express them in 
terms of the 2-dimensional connection A % a and the triad field Ej — e bc e^i]ji, where a = 1,2 are space coordinate indices 
and i, j = 1, 2, 3 are su(2) indices. The symplectic structure is defined by 

{A i a {x),E){y)}=5 b a 5 i j 8^{x,y). 

Local symmetries of the theory are generated by the first class constraints D^E^ ~ and F* b (A) ~ 0. 

To quantize the theory we first find a representation of the basic variables in an auxiliary Hilbert space 7i a ux- The 
basic functionals of the connection are represented by the set of holonomies along paths 7 C £. Given a connection 
A and a path 7, one defines the holonomy ft, 7 [yl] by 

hy[A] = Fcxp / A. (81) 

J -y 

As for the triad, its associated basic variable is given by the smearing of E along co-dimension 1 surfaces. One 
promotes these basic variables to operators acting on an auxiliary Hilbert space where constraints are represented. 
The physical Hilbert space is defined by those 'states' that are annihilated by the constraints. As these 'states' are 
not normalizable with respect to the auxiliary inner product they are not in H aux and have to be regarded rather as 
distributional. 

The auxiliary Hilbert space is defined by the Cauchy completion of the space of cylindrical functionals Cyl, on the 
space of (generalized) connections A. The space Cyl is defined as follows: any element of Cyl, ^rj[A] is a functional 
of A labeled by a finite graph T C S and a continuous function / : SU(2) Ne ^ — > C where A^(r) is the number of 
links of the graph T. Such a functional is defined as follows 

9 r ,f[A]=f(h n [A],---,l hst(n [A]) (82) 

where h li [A] is the holonomy along the link 7, of the graph T. If one considers a new graph F' such that V C V, then 
any cylindrical function $rj[A] trivially corresponds to a cylindrical function ^r',/'^]- This space can be endowed 
with the well known Ashtekar-Lewandowski inner product. The Gauss constraint can be defined in terms of the 
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basic variables introduced above. It generates gauge transformations whose action on Cyl transforms the holonomy 
as follows 

f h [A]^g.? h [A]gt 1 (83) 

where g s ,gt & SU{2) are group elements associated to the source and target nodes of 7 respectively. The so-called 
kinematical Hilbert space Hkin C H a ux is defined by the states in 7i aux which are gauge invariant and hence in the 
kernel of the Gauss constraint. 

A basis of gauge invariant functions can then be constructed by contracting the tensor product of representation 
matrices with sw(2)-invariant tensors or su(2)-intertwiners. If we select an orthonormal basis of intertwiners t„ £ 
Inv[ji (g> J2 <8> ■ ■ • <8> jwj where n labels the elements of the basis we can write a basis of gauge invariant elements of 
Cyl called the spin network basis. Each spin network is labeled by a graph r C S, a set of spins j 7 labeling links 7 
of the graph T and a set of intertwiners i n labeling nodes n of the graph T, namely: 

Sr, iiehM [A] = ® in ® S (K[A]) • (84) 

ner 7er 

In order to lighten notations, we will omit indices (the graph, representations and intertwiners) labeling spin-networks 
in the sequel. One can take the kinematical inner product of two spin networks s, s', which we will denote < s\s' > 
by integrating expression (|84p with the Ashtekar-Lewandowski measure. 

Let us review the group averaging procedure to generate the projector, following closely the discussion of Noui and 
Perez. One starts with the formal expression 

P = J] 5(F(A))» = J D[N] exp(i / Tr[NF(A)]) , (85) 

where N £ su(2). We now introduce a regularization of ([55]) . We will give a definition of P by providing a regular- 
ization of its matrix elements < Ps, s' > for any pair of spin network states s,s' S Ti-Mn- Let's denote by T and V 
the graphs on which s and s' are defined respectively. We introduce an arbitrary cellular decomposition of S denoted 
r , where eg J?, such that the graphs T and V are both contained in the graph defined by the union of 0-cells and 
1-cells in r and for each individual 2-cells (plaquette) p there exist a ball B € of radius e — defined using the local 
topology — such that pc6 ( . 

Consequently all 2-cells shrink to zero when e — > 0. We consider a local patch [/cS where we choose the cellular 
decomposition to be square with cells of coordinate length e. In that patch, the integral in the exponential in 
can be written as a Riemann sum 



F[N] = J Tr[NF(A)] = lim ^ e 2 Tr[A^F pI ], 



p 1 

:th 



where p % labels the i th plaquette and N p i £ su(2) and F p i £ su(2) are values of N^Tj and Tje ab F^ b [A] at some interior 
point of the plaquette p 1 and e ab is the Levi-Civita tensor. The basic observation is that the holonomy U p i £ SU (2) 
around the plaquette p % can be written as 



which implies 



U pi [A] = l + e 2 F p ,(A)+0{e 2 ) 
F[N] = J Ti[NF(A)} = lim ^Tr[JVpiJ7 pl [A]] , (86) 



where the Tr in the right hand side is taken in the fundamental representation. Notice that the explicit dependence 
on the regulator e has dropped out of the sum on the r.h.s., a sign that we should be able to remove the regulator 
upon quantization. With all this it is natural to write the following formal expression for the generalized projection 
operator, 

<Ps,s'> = lim <Il/ dN P' exp(iTr[A^l> pI ])s, s' > (87) 

p i 

= lim <Y[ S (Up*)s, s' >, 
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where the last equality follows from direct integration over N p t at the classical level; 8(U) being the distribution such 
that J dg f{g)6(g) — /(l) for / e C 2 (SU{2)). We can write 8{U) as a sum over unitary irreducible representations of 
SU(2), namely 8{U p i) = + 1) Xj(U p i), where Xj{U) is the trace of the j-representation matrix of U € SU(2). 

It is therefore straightforward to evaluate explicitly the action of the projector. 

We now proceed to construct the projector using our technique. We construct a Hamiltonian 



where F p T l is the element of the SU (2) algebra such that exponentiated gives the holonomy around the plaquette p 
(t 1 are the Pauli matrices), U(p) = exp(iF p ). We now proceed to construct the projector as usual, 

V= lun ( (89) 

A^oo V 7T / 

If this expression is to agree with the projector computed by Noui and Perez, one has to have that, 

x^L J (n d9 ^J f(9pi,9 P 2,---)'Px = J (jj^dgij f{9p 1 ,9p 2 ,---)'[[8(g p ). (90) 



In order to prove this identity we start by introducing a parameterization of the SU(2) group [2l| . 

a b 
-b* a* 



where 



b = = — i sin 



(|) -i sin (|) cos 0, (92) 
(|)sin^. (93) 



One then has that 



and 



J dg = 4 J sin 2 (^j dw J sinOdd J dip, 



(94) 



X» = V 7-jt ' ■ (95) 



And we see that u> = corresponds to a = 1, b = 0. The -F's take the form, 



From here we can compute H, 



and the projector can be written as, 



F 1 = - ^-J sin cos <p, (96) 

F 2 = - (^j sin0sin</>, (97) 

F 3 — ~ cos0. (98) 

H = J2<4, (99) 



■ \ \ 3/2 

iX x 



"-fitn^j e """ ; ' <100) 
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To complete the proof we would like to show that 



s (9 P ) = ^^s(u) = s(g P (u,o,<p)), (ioi) 



which can be shown by noticing that 

dg P f{9p)5(g p {uJ,9,ip)) = 2 J du> J dftsin 2 (~J f(g p {u,6,ip)) — JjTj^H 



2 J* dcsin 2 (f) f<M P (wAv>))-£sS(w) 
2 dujf (g p (uj, 6, V))<SH^ sin 2 (|) 

doif(g p (u, 9, <p))S(u) cos - /(l). (102) 



Finally, using the following Dirac identity for distributions, 



X --E(g) ^K) 



e *=2J^J 77T-Vf7^K)> (103) 



where <5^ 2 ^ is the (2j)-th derivative of the Dirac delta. Inserting this expression in poop and using (|102[) we prove 
the identity (|90[) which implied that the projector we constructed with our technique agrees with the one obtained 
via group averaging by Noui and Perez. 



VIII. CONCLUSIONS 



We have introduced a new approach to the quantization of totally constrained systems and worked out several 
examples. One reaches two conclusions. On one hand, one can use the technique to construct in a well defined way 
the physical space of states of the totally constrained theory of interest. Of special appeal is the fact that one does 
it without the need of promoting the constraints to quantum operators, which in several situations is known to be 
problematic. In particular one does not have to deal with the problem of reproducing the quantum constraint algebra. 
In examples we have shown that the results coincide with those of the group averaging procedure, where the latter 
exists. Moreover one has a correspondence principle in that one needs to reproduce the classical equations of motion 
as quantum operatorial relations that offers guidance in the intermediate steps of the process. The technique therefore 
has aspects that arc quite distinct from conventional ones, especially when one treats field theories. 

On the other hand the technique can be viewed as a new paradigm for dealing with cases where the continuum 
theory does not exist. In those cases one builds discrete theories where one can introduce a relational notion of 
time and compute probabilities and expectation values. In examples we have shown that if one is able to take the 
continuum limit, one reproduces predictions of the evolving constants approach and therefore recovers all the physics 
in the theory. But we have also shown in examples where the continuum limit does not exist that one can define 
semiclassical regimes that reproduce, at certain scales, the classical results of the theory of interest. Although we do 
not yet know if the continuum limit exists in the case of general relativity, this point of view would agree with popular 
beliefs in that case: one would have a fundamental, discrete theory that reproduces, at large scales, general relativity 
even if the quantum theory strictly does not exist in the continuum limit. 

Finally, the uniform discretizations also limit importantly the usual ambiguities that appear when one discretizes 
theories and may even be useful in the classical modeling of constrained theories. Future steps will require to test the 
technique in cases of increasing complexity in the gravitational case. We have shown that in the cosmological case 
one reproduces the results of loop quantum cosmology and the technique may allow to explore some of the attractive 
features of the latter in more complex models in the near future. 
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